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The cosmological dynamics of minimally coupled scalar field that couples to the back- 
ground matter with thermal interactions is investigated by using statefinder diagnostics. 
The time evolution of the statefinder pairs {r, s} and {r, q} are obtained under the 
circumstance that different values of model parameters are chosen. We show that the 
thermal coupling term does not affect the location of the late-time attractor, but exert 
an influence on the evolution of the statefinder parameters. The most notable feature of 
the r — s plane for the thermal coupling model which is distinguished from the other dark 
energy models is that some part of the curve with thermal coupling can form a closed 
loop in the second quadrant (r > 1, s < 0). 
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1. Introduction 

Various of cosmological observations strongly suggest that our universe is undergo- 
ing an accelerated expansion phase. In order to explain the acceleration, an unex- 
pected energy component of the cosmic budget, which dominate the universe only 
recently, is introduced by many cosmologists. However, the physical origin of dark 
energy is still mysterious. 

Perhaps the simplest proposal is the Einstein's cosmological constant A (vac- 
uum energy), whose energy density remains constant with time. However, due to 
so me con ceptual problems associated with the cosmological constant (for a review, 
see HEEl) ; a i ar g e variety of alternative possibilities have been explored. The most 
popular among them is quintessence scenario which uses a scalar field <fi with a 
suitably chosen potential V((j>) so as to make the vacuum energy vary with time. A 
form of quintessence called " tracker fields" , whose evolution is largely insensitive to 
initial conditions and at late times begin to dominate the universe with a negative 
equation of state, was introduced to avoids the problems relevant to cosmological 
constant SI. Another approach to solve the puzzle is introducing an interaction term 



f 
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in the equations of motion, which describes the energy flow between the dark energy 
and the rest matter (mainly the dark matter) in the universe. It is found that, with 
the help of a suitable coupling, it is possible to reproduce any scaling solutions. In- 
clusion of a non-minimal coupling to gravity in quintessence models together with 
further generalization leads to models of dark energy in a scalar-tensor theory of 
gravity. Besides, some other models invoke unusual material in the universe such as 
Chaplygin gas, tachyon, phantom or k-essence (see, for a review, ^ and reference 
therein) . 

As so many dark energy models have been proposed, it becomes urgent to give 
them an unambiguous discrimination. The equation of state w, could discriminate 
some basic dark energy models properly, for example, w = — 1 for cosmological 
constant, w > — 1 for quintessence and w < — 1 for phantom, but as more and more 
interacting models were investigated, the ambit of w become not so clear. A new 
geometrical diagnostic, dubbed the statefinder pair {r, s} is proposed by Sahni et 
al where r is only determined by the scalar factor a and its derivatives with 
respect to the cosmic time t, just as the Hubble parameter H and the deceleration 
parameter q, and s is a simple combination of r and q. The statefin der pair has 
been used to explor a series of dark energy and cosmological models I 7 | 8 | 9 | 10 | As 
is analyzed, the "distance" from a given dark energy model to the LCDM scenario 
can be clearly identified via the r(s) evolution diagram. The current values of r 
and s, which can be calculated in models, are evidently valuable since it is expected 
that they can be extracted from data coming from SNAP (SuperNovae Acceleration 
Probe) type experiments. Therefore, the statefinder diagnostic combined with future 
observ ation s may possibly be used to discriminate between different dark energy 



The effects of thermal coupling between the quintessence field and the ordi- 
nary matter particles was first investigated by Hsu and Murray H| They set the 
quintessence field to be a static external source for a Euclidean path integral de- 
picting the thermal degree of freedom and let the time-like boundary conditions of 
the path integral have a period which is decided by the inverse of temperature. And 
they show that i f in the early universe matter particles are in thermal equilibrium, 
quantum gravity will induce an effective thermal mass term for quintessence field 
<fi, which takes the form 



where Ta o _1 denotes a measurement of temperature of background matter, Mp 
the Planck mass scale, (3 a dimensionless constant, a the scale factor (with current 
value a$ = 1), and find that even Planck-suppressed interactions between matter 
and the quintessence field can alter its evolution qualitatively. In a previous paper 
HSI, we investigated the dynamics of the cosmology with quintessence using the 
phase space analysis, which has the above thermal coupling to the matter in a 
complete manner and analyzed the conditions for the existence and stability of 



models^. 




(1) 
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various critical points as well as their cosmological implications. In this paper, we 
apply the statefinder diagnostic to the cosmological dynamics of minimally coupled 
scalar field that couples to the background matter with thermal interactions. 

2. The Setup 

Let us study quintessence with thermal coupling to ordinary matter particles in 
spatially flat FRW cosmological background 

ds 2 = dt 2 - a 2 {t)dx 2 . (2) 

For the spatially homogeneous scalar field minimally coupled to gravity with thermal 
interaction ([J), the evolution is governed by the Klein-Gordon equation 

+ 3i^ + y'(^)=2(-^ ^t 4 , (3) 

where the overdots denote the derivative with respect to cosmic time and the prime 
denotes the derivative with respect to (f>. Here the Hubble parameter H = a/a is 
determined by the Friedmann equation 



and 



Pm + \ft + V((f>) 



(4) 



H = - — {pm + P,n + 4>), (5) 

where k 2 = 8ir/M 2 , p rn and p m are the energy density and pressure of the bary- 
otropic matter, respectively. From Eqs.([3])-([5]) and the conservation of energy, p m 
satisfies the following continuous equation 

p m + 3H(p m +p m ) = -2(-^) WT\ (6) 

and p m = (7— l)/3m, where 7 is a constant, < 7 < 2, such as radiation(7 = 4/3) or 
dust (7 = 1). It is clear that when the thermal coupling parameter j3 becomes zero, 
the equations (J3j- ((6j will return to those of the standard one scalar field quintessence 
scenario. In current situation, both quintessence and baryotropic matter are not 
conserved, but the conservation of the overall energy holds. 

To be concrete, in this paper we only consider the quintessence with an expo- 
nential potential energy density, i.e., 

V{<t>) = V exp(-X K<P) (7) 

where the parameter Vq and Ao are two positive constants. Exponential potentials 
have been studied extensively in various situations, and these are of interest for two 
main reasons. Firstly, they can be derived from a good candidate of fundamental 
theory for such being string/M theory; secondly, the equations of motion can be 
written as an autonomous system in the situation. 
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Introducing the following dimensionless variables: 



K (j) nyjv{4>) k 4> V6 

' y= at tt ' 2 = 77 — aT7T7> £ = _ T> ^ = lna (8) 



y/6H' ' M P V3a 2 7j' k0 

and taking the potential |J7J| into account, we can rewrite the equation system (J4jl- (J6J) 
as the following autonomous system: 

% = |x[ 7 (l -x 2 - y 2 ) + 2a; 2 ] - (3. + /3V£ - ±A 0j/ 2 ), 
= ^[7(1 - z 2 - y 2 ) + 2a: 2 ] -2z + xz^, 

3. Statefinder Parameters 

The traditional geometrical diagnostics, i.e., the Hubble parameter H and the de- 
celeration parameter q = —da/ a 2 , are two good choices to describe the expansion 
state of our universe but they can not characterize the cosmological models uniquely, 
because a quite number of models may just correspond to the same current value 
of H and q. Fortunately, as is shown in many literatures, the statefinder pair {r, s} 
which is also a geometrical diagnostic, is able to distinguish a series of cosmological 
models successfully. 

The statefinder pair {r, s} defines two new cosmological parameters in addition 
to H and q: 

1 d 3 a r — 1 

r= ~alP~aW , S= 3(g-l/2)' ( ' 

As an important function, the statefinder can allow us to differentiate between 
a given dark energy model and the simplest of all models, i.e., the cosmological 
constant A. For the ACDM model, the statefinder diagnostic pair {r, s} takes the 
constant value {1,0}, and for the SCDM model, {1, 1}. 

From Eqs. (fTU)) and ©, the statefinder parameter r can be explicitly written as 

r = 1 + ^7(7 - 1)(1 - x 2 - y 2 ) + 3(7 - 2)0 2 xz 2 £ > + 9x 2 - 3X xy. (11) 
And the deceleration parameter q is given by 

q=~7(l-x 2 -y 2 )+3x 2 -l. (12) 

Therefore, the other statefinder parameter s is equal to 

2 -y 2 ) + f(7-2)/3 2 xz 2 £ 
7(1 - x 2 - y 2 ) + 2x 2 - 1 



7(7 - 1)(1 - x 2 - y 2 ) + f( 7 - 2)fP-xz 2 i + 2x 2 - \\ xy 



(13) 
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We draw a 7 — A 2 plane in Figure Q] and divide it into different regions. From 
the analysis of our previous paper we know that the equation system, Eqs.©, 
have a sole late-time attractor only if the conditions 7 < 4/3 and/or A 2 < 24 are 
satisfied. It is easy to find that these conditions correspond to region A and B in 
Figure [T] and the common boundary of region A and B is the segment of straight line 
Aq = I87. In what follows, we call the condition that region A satisfies case A and 
that region B satisfies case B. For case A and case B, the attractor (x c ,y c , z Cl ^ c ) 
is located at the point (Ao/6, ^/36 — Aq/6, 0, 0) and (37/A0, 3-^/7(2 — 7)/ Ao, 0, 0) in 
the phase space, respectively. It is readily to check that in the common boundary 
of region A and B, there actually exists only one attractor as is expected. It is also 
worth pointing out that the attractor depends only on the parameter A for case A, 
and for both cases, the attractor has nothing to do with the coupling parameter (3. 
In statefinder parameter space, the attractor lies at the point (r c , q c , s c ) and it can 
be expressed explicitly by 



A, 



Tc = 1+ 12 ' x ' m " A 

\2 

9c = -l + T^ 



12 

3 - ^36 - Aq 



for case A and 



= 1 + 777(7-1) + 



3 (A 2 -18) 

81(7 - 2) 7 2 277v/(2-7)7 



(14) 



2 ' v ' ' A 2 A 



3 7 

- 1 + T< 



18(7 - 2) 7 2 6 7v /(2-7)7 
(1-7)A§ (l-7)Ao 



7+— ZXT2- + — n IT - ( 15 ) 



for case B, respectively. 

Let us now turn to statefinder parameters plane analysis. We shall investigate 
the properties of the evolution trajectories of statefinder parameters with three 
parameters (3, 7 and Ao, by fixing two of them and varying the rest. Note that 
we are only interested in case A and case B, due to the existence of the attractor 
solution. 

We first fix the brayotropic parameter 7 and potential parameter Ao, and let 
the thermal coupling constant /3 be adjustable. We show the time evolution of the 
statefinder pairs {r, s} (left panels) and {r, q} (right panels) in Figure [21 Figure 
[3] and Figure |4] where the fixed parameters 7 = 1 and Ao = 1,4,7, respectively. 
In these three figures, we choose an initial state that the universe was ongoing a 
decreasing expansion and the kinetic energy of the quintessence is negligible. The 
trajectories plotted in Figure [2] and Figure [3] correspond to case A. Obviously, they 
are all located in the second (s < 0, r > 1) and/or fourth (s > 0, r < 1) quadrants in 
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2/3 i 4/3 7 

Fig. 1. The schematic illustration of the parameter space. In region A and B, the autonomous 
system J9j has a late-time attractor, respectively. 



the r — s plane. This is not surprising because we have assumed that the equation of 
state of quintessence field became negative after some past moment and this makes 
the deceleration parameter q < 1/2 from then on as is shown in right panels of 
Figure [5] and Figure [31 From the definition of statefinder parameter s, Eq. lfTU)) . it 
is easy to find that if r < 1 then s > 0, if r > 1 then s < and if r = 1 and 
q ^ 1/2 then s = which is just corresponding to the ACDM fixed point. It is 
worth pointing out that the attractors in the two figures have different cosmological 
implication: the attractor illustrated in Figure [2] means eternal acceleration of the 
universe, while the one in Figure [3] means the universe will decelerate in the future. 
The trajectories plotted in Figure [4] correspond to case B. The accelerating rate of 
expansion of the universe will be determined only by the baryotropic matter. We 
choose the baryotropic parameter 7 — 1 in Figure [4] which means the universe will 
be decelerated in the future, as is shown in the r — q plane (right panel). From 
these three figures, we see that the thermal coupling term enhances the possible 
maximum value of r and makes it first increase to a maximum value (> 1) then 
decrease. This effect reflected in the r — s plane is that some parts of the (3 ^ 
curves can be located in the second quadrants (r > 1 , s < 0) and these parts form 
closed loops. This feature of the r — s plane for the thermal interacting quintessence 
model is distinguished from those of ordinary quintessence and other dark energy 
models. 

Next, let us show the dependence of the time evolution of the statefinder pair 
{r, s} on the brayotropic parameter 7 and potential parameter Ao in the absence or 
presence of the thermal interaction. In Figurc[5l without loss of generality the bray- 
otropic parameter 7 is fixed to be unity and the potential parameter is selected to 
satisfy the condition Ag = 12, 18, 24, respectively. The left panel shows the evolution 
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Fig. 2. The time evolution of the statefinder pairs {r, s} and {r, q} is illustrated in left and right 
panel, respectively. In both panels the model parameter 7 and Ao are fixed to be unity and the 
thermal coupling constant /3 takes different values. For (3 = curve, r and s monotonically decrease 
from 1 and increase from to the value the attractor takes, respectively, and q also monotonically 
decrease from positive to negative. For j3 ^ curves, q keeps monotonically decreasing, whereas r 
(s) first increase (decrease) from 1 (0) to a maximum (minimum) value then decrease (increase) 
and pass through the ACDM fixed point (in r — s plane) to the attractor. 




-4-3-2-101234 " -3.6 -0.4 -0.2 0.2 0.4 



Fig. 3. The time evolution of the statefinder pairs {r, s} and {r, q} is illustrated in left and right 
panel, respectively. Here the fixed model parameters are chosen to be 7 = 1 and Ao = 4, whereas 
the thermal coupling constant takes different values. For /3 = curve, s monotonically increase 
from to the value the attractor takes, while r and q first decrease from the start point to a 
minimum value then increase to the attractor. For f3 ^ curves, q has the similar behavior of that 
of /3 = curve, whereas r (s) first increase (decrease) from 1 (0) to a maximum (minimum) value 
then decrease (increase) and after passing through the ACDM fixed point (in r — s plane), they 
show a similar behavior of that of ft = curve. 

for the case that thermal interaction is absent (/3 = 0) and the right one for the case 
with thermal interaction (we set (3 = 4 here for example). Clearly, the effect that 
the thermal coupling term exerts is the same as that pointed out previously. The 
Aq = 12 curve has a definite end point which locates the late-time attractor but the 
others haven't. This is not surprising. The parameter related to the Aq = 12 curve 
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Fig. 4. The time evolution of the statefinder pairs {r, s} and {r, q} is illustrated in left and right 
panel, respectively. Here the fixed model parameters are chosen to be 7 = 1 and Ao = 7, while the 
thermal coupling constant j3 takes different values. Because q approaches to and takes the value 
of 1/2 many times at the late time, but at the same time r is always less than unity, the value of 
s will become infinity many times. Therefore, the plot in the left panel shows only a piece of the 
whole time evolution. 

lies in the region A of the parameter space (see Figure [1]) and the location of the at- 
tractor is determined by Eqs. (fT4|) which take definite values, while the parameters 
relevant to the other curves locate in the region B or the boundary line between 
region A and B, so the locations of the corresponding attractor is determined by 
Eqs. (|15p which are infinite for the parameters selected. The situation is similar in 
Figure where we fix the value of A and select three different representative values 
for 7. In addition, the dark dots marked on the curves in the two figures represent 
the present values of the statefinder parameters {r , so}- It should be noted that the 
true values of {tq, sq} of the universe should be determined in a model- independent 
way, which we can only place our hope on the future experiments. In principle, if 
this is achieved, the statefinder parameters can be viewed as a discriminator for 
testing various cosmological models. 

4. Conclusions 

We have investigated the time evolution of quintessence field with an exponential 
potential and with or without thermal coupling to the background matter in a 
flat FRW universe from the statefinder viewpoint in this paper. The statefinder 
diagnostic have been performed. It is found that that the thermal coupling term 
does not affect the location of the late-time attractor which is determined only by 
the potential of the quintessence field and the background baryotropic matter, but 
exert an influence on the evolution of the statefinder parameters. The most notable 
feature for the thermal coupling model which is distinguished from the other dark 
energy models is that the curve with thermal coupling can form a closed loop in 
the second quadrant (r > 1, s < 0) in the r — s plane. What we expect now is that 
more and more accurate observational data could be offered to determine the model 
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Fig. 5. The trajectories in the statefinder plane {r, s} for the model where we fix the baryotropic 
parameter 7 to be unity and let the thermal coupling constant (3 = and /3 = 4 in the left and 
right panel, respectively. The dark dots indicate the current value of statefinder parameters. The 
triangle at the end of the Ag = 12 curve denotes the location of the late-time attractor. For the 
other two curves the attractor is located at infinity. 




Fig. 6. The trajectories in the statefinder plane {r, s} for the model where we fix the baryotropic 
parameter Ao = %/l8 and let the thermal coupling constant /3 = and /3 = 4 in the left and 
right panel, respectively. The dark dots indicate the current value of statefinder parameters. The 
triangle at the end of the 7 = 2/3 curve denotes the location of the late-time attractor. For the 
other two curves the attractor is located at infinity. 



parameters more and more precisely, rule out some models and consequently shed 
light on the essence of dark energy. 
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